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OPERATOR SPLITTINGS AND SPATIAL APPROXIMATIONS FOR 

EVOLUTION EQUATIONS 
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To Rainer Nagel, our teacher and friend. 


Abstract. The convergence of various operator splitting procedures, such as the se¬ 
quential, the Strang and the weighted splitting, is investigated in the presence of a 
spatial approximation. To this end the relevant notions and results of numerical analysis 
are presented, a variant of Chernoff’s product formula is proved and the general Trotter- 
Kato approximation theorem is used. The methods are applied to an abstract partial 
delay differential equation. 


1. Introduction 

Operator splitting procedures are usually used to solve partial differential equations 
numerically. They can be considered as certain time-discretization methods which simplify 
or even make the numerical treatment of differential equations possible. 

The idea behind operator splitting procedures is the following. Usually, a certain physical 
phenomenon is the combined effect of several processes. The behaviour of a physical 
quantity is described by a partial differential equation in which the local time derivative 
depends on the sum of the sub-operators corresponding to the different processes. These 
sub-operators usually are of different nature. For each sub-problem corresponding to each 
sub-operator there might be an effective numerical method providing fast and accurate 
solutions. For the sum of these sub-operators, however, we usually cannot find an adequate 
method. Hence, application of operator splitting procedures means that instead of the 
sum we treat the sub-operators separately. The solution of the original problem is then 
obtained from the numerical solutions of the sub-problems. A simple splitting procedure 
was proposed by Bagrinovskii and Godunov (see [2|) in 1957 as an example. However, they 
were systematically studied only in 1968 by Marchuk (see j20[) and Strang (see (25[,j26|). 
Since then operator splitting procedures have been widely applied to various physical 
processes, see e.g., Zlatev [28=], Zlatev and Dimov [29=], or Botchev et al. [3|. 

Although operator splitting combined with spatial approximations is widely used in appli¬ 
cations, the theoretical convergence analysis of these problems seems to be still missing. 
One aim of this paper is to investigate systematically the convergence of various splitting 
procedures combined with spatial approximations. This important question is motivated 
by numerical analysis and real-life applications but in the present paper it will be inves¬ 
tigated in the framework of functional analysis. This means that we consider and apply 
the well-known theorems of Chernoff and Trotter-Kato from another point of view. They 
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are just mathematical tools to prove the convergence of numerical schemes solving partial 
differential equations. 

First, we recall various notions from numerical analysis concentrating on finite difference 
methods and, in particular, the sequential, Strang, and weighted splitting is introduced. 
In Section 3 we give an overview on the convergence of these splitting procedures, and the 
convergence is investigated in the presence of a spatial approximation. The last section is 
devoted to the application of these results to delay equations. Further applications along 
with numerical experiments and error estimates will be presented in a subsequent paper. 
Since we apply the results of operator semigroup theory, our main reference regarding 
notations and terminology is Engel and Nagel ( 9 ]. To avoid technical complications, we 
restrict ourselves to the case of two sub-operators. The general problem can be treated 
analogously. Further, we make the following assumption for the rest of this work. 

Assumption 1.1. Assume that A is a Banach space, A and B are closed, densely defined 
linear operators generating the strongly continuous operator semigroups T = (T(i)) 

and S = (S(t)) t>0 , respectively. Further, we assume that the closure of A + B, A + B 
with D(A) (~) D(B) C D(A + B) is also the generator of a strongly continuous semigroup 

u =kww 


2. Finite difference methods and splitting procedures 

In this section we recall some basic notions of numerical analysis needed for our later 
investigations. Though the results presented are standard and well-known, it seemed to 
be appropriate to present them to a reader having his interests in evolution equations 
and operator theory. Since we use the finite difference method for solving the equations 
numerically, we restrict ourselves to some of its important properties. Our discussion 
follows Section 6 of Atkinson and Han |lj, but another relevant reference is Richtmyer 
and Morton j24j. Then we introduce three types of splitting procedures investigated in 
the present work: sequential, Stang, and weighted splittings. We will see that splitting 
procedures are special finite difference methods. 

2.1. Finite difference methods. Finite difference methods approximate the differential 
operator by a difference operator. Thus, instead of the differential equation, a system 
of algebraic equations has to be solved at each time step. A numerical scheme is only 
applicable if it can provide numerical solution arbitrarily close to the exact solution, i.e., 
if the method is convergent. In what follows, we define the basic concepts of numerical 
analysis we need later on. 

Let X be a Banach space and G : D{G) C X —y X be a linear (usually unbounded), closed, 
densely defined operator. Consider the abstract Cauchy problem, also called initial value 
problem 

d u(t) 

( 1 ) 


= Gu(t), t > 0, 


dt 

w(0) = x G X. 

Recall that a function u : M + —y X is called a (classical) solution of the abstract Cauchy 
problem (0Q) if the function u is continuously differentiable, u(t) G D(G), and 

u{t + h) — u(t) 


lim 

h-> o 


— Gu(t ) 


= 0 
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with u(0) = x. This limit is understood to be the right limit at t — 0. Since G generates 
the strongly continuous semigroup (i7(i)) , it is clear that problem (JTJ) has a unique 

solution of the form 

(2) u(t) = U(t)x for all t > 0, x G D{G). 

Although the solutions of (JTJ) are defined for all t > 0, with the help of a computer we can 
only hope for a solution in finite time. Therefore, from now on we choose an arbitrary but 
fixed “end point” t 0 G R + , that is, we solve ([I]) numerically on the time interval [0, £ 0 ]. 

Definition 2.1. A finite difference method is a one-parameter family of linear operators 

F(h) : X — > X, where h G (0,t 0 ], f 0 G R + fixed. 

The approximate solution u m (t ) at time t = mh is then defined by 

u m (t) = F{h) m x = F(t/m) m x for m G N. 

We usually refer to h as the time step of the numerical method. 

Now we can define the consistency, the convergence, and the stability of a difference 
method. 

Definition 2.2 (Consistency). The difference method is called consistent if for all x G 
D{G) and for the corresponding solutions u of the abstract Cauchy problem (JT]) the 
following holds: 

F(h)u(t) — u(t + h) _ 
h 

uniformly for t G [0,to]. 

Roughly speaking, consistency means that the approximating difference equations con¬ 
verge in some sense to the original abstract Cauchy problem. 

Definition 2.3 (Convergence). The difference method is called convergent at time t G 
[0, to] if f° r fixed t and for any x G X we have 

(3) lim \\F(hi) mi x-u(t)\\ = 0, 

hi~>0 

where (m*)^ is a sequence of integers and is a null-sequence of step sizes such 

that mihi = t. 

Remark 2.4. Convergence means that the numerical solution tends to the exact solution 
of the problem as the time step tends to zero. We defined the convergence only at a fixed 
time level. However, from the theorems presented later it follows that the convergence is 
uniform for t in compact intervals. 

Definition 2.5 (Stability). The difference method is called stable if the family of opera¬ 
tors 

{F(h) m : h G (0, i 0 ], m h < t 0 , m G N} 
is uniformly bounded, i.e., there exists a constant M > 0 such that 

||F(h) m || < M for mh < to and for all h G (0,to]. 


lim 

o 


In other words, stability means that the approximating solutions remain bounded. 
The above three notions are connected through the following result (see Lax 
8 in Section 34.3] and Atkinson and Han jl;, Theorem. 6.2.11]). 


19, Theorem. 
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Theorem 2.6 (Lax Equivalence Theorem). Let us assume that the initial value problem 
m is well-posed. Then a consistent finite difference scheme is convergent if and only if it 
is stable. 


2.2. Splitting procedures. In order to introduce the operator splitting procedures, we 
consider the following abstract Cauchy problem on the Banach space X : 


(ACP) 


^ = (A + B)u(t), t> 0, 
w(0) = x £ X, 


where the operators A and B satisfy Assumption 11.11 

Since operator splitting procedures are time-discretization methods (see the explanation 
below), analogously to the numerical time step introduced in Definition 12.11 we choose 
the splitting time step of length r £ M + . 


Definition 2.7. The solution obtained by applying a splitting procedure is called split 
solution. We remark that the split solution is only defined on the mesh 

(4) u; r := {kr, k £ M}. 


In the following we collect three splitting, procedures (see, e.g., Bagrinovskii and Godunov 

and consider the corresponding sub-problems 


Q, Strang [25[, Marchuk 20|, Farago 
with their solutions. Recall that we assumed that (A,D(A)) and (5, D(5)) generate 
the strongly continuous semigroups (T(i)) and (S'(t)) , respectively (cf. Assumption 

in]). 


2.2.1. Sequential splitting. As we have already mentioned in the Introduction, applying 
splitting procedures means that we solve sub-problems related to the sub-operators A and 
B separately. For sequential splitting this process can be formulated as follows. 

| = fe((fe-l)r,fcr], 

[ u[ k \(k - 1)t) = u sq ((k - 1)t), 

fe((A-l)r,fer], 

u { 2 k \(k - l)r) = u[ k \kr), 
u sq (kr) := u^^kr), 

with f 6 N and u sq (0) = x. Using that u[ k \t) = T(t)u sq ((k — 1 )r), and that u^\t) = 
S(t)u ( {\kT ) = S(t)T(r)u sq ((k — l)r), we see by a simple induction argument that the 
split solution u sq (kr), obtained by applying the sequential splitting procedure, can be 
written as 

(5) u sq (kr ) = [5’(r)T(r)] fe x for k £ N, kr < t 0 , and x £ X. 

Since this is a finite difference method, the convergence of a splitting procedure also 
needs to be investigated at a certain time level. Therefore, we write formula (J5]) in a more 
convenient way. Now we do not fix the value of the splitting time step r, but fix the value 
of t > 0. With r \—tjn we obtain the solution 

(6) < q (f) = [S(t/n)T(t/n)} n x 
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for all n G N, t > 0 and x G X. Hence, sequential splitting is a finite difference method 
with 

F sq (h) = S(h)T(h), h G (0, t 0 ]- 

2.2.2. Strang splitting. In the case of this splitting technique three sub-problems have to 
be solved for one splitting time step: 

| dU | 1 ^ ) = Au[ k) {t), t G ((k - l)r, (k - \) t\ , 

[^((/c-ljr) = u St ((k — 1 )t), 

| = Bu}\t), t G ((k - 1)t, kr], 

\ u}\(k - 1)t) = u[ k) ((k - \) t) , 

—= Auf\t), te ((k — |) r, fcr] , 

4 fe) (( fc - |) T ) = ^(kr), 
u s \kr) :=u {k \kr), 

where u St (0) = x and k G N. The split solution can be written in this case, using the 
same arguments as above, as 

(7) u s \kr) = [T(t/2)S'(t)T(t/2)] A: x for k e N and x E X. 

As above, substituting r := t/n with t > 0 hxed, we have 

(8) nf(i) = [T(t/2n)S(t/n)T(t/2n)] n x 

for all n G N, t > 0 and x e X. Hence, Strang splitting is a finite difference method with 

F st (h ) = T(h/2)S(h)T(h/2), h G (0, t 0 ]. 

2.2.3. Weighted splitting. It is obtained by using two sequential splittings: first starting 
with operator A, and then starting with operator B. At time t — kr the split solution 
is computed as a weighted average of the split solutions obtained by the two sequential 
splitting steps: 

u w (kr) = Qu sq,AB (kr) + (1 - Q)u sq ’ BA (kr), 

where 0 G (0,1) is a given weight parameter, and u sq,AB [kr ) and u sq ’ BA (kr ) are the split 
solutions of the above two sequential splittings at time kr, respectively. In this case the 
split solution has the form 

(9) u w (kr ) = [&S(t)T(t) + (1 - 0)T(r)^(r)] fc x 

for k G N and x G X. With varying splitting time step r := t/n and fixed t > 0, the 
above formula can be rewritten as 

(10) = [ QS(t/n)T(t/n ) + (1 — Q)T(t/n)S(t/n)] n x 

for all n G N, t > 0, and x G X. Hence, weighted splitting is a finite difference method 
with 

F w (h ) = QS{h)T{h) + (1 - G)T{h)S{h), h G (0, t 0 ]. 

The case 0 = \ is called symmetrically weighted splitting. 
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Remark 2.8. Let us make some notes on the various splitting procedures. The consistency 
(and its order) of the different splittings were shown by Bjprhus in [h| and Farago and 
Havasi in [l3[ for general Co-semigroups under a restrictive domain condition on the 
generators. For matrices it is well-known that the sequential splitting is of first order, the 
Strang and the symmetrically weighted splitting is of second order. The weighted splittigs 
were introduced because they allow the use of parallel computations. For analytic and 
unitary (semi)groups in Hilbert spaces, the convergence of sequential and Strang splittings 
were also investigated see, e.g., Ichinose et al. 14||, Neidhardt and Zagrebnov [2l[ and [22j |. 


and Zagrebnov [27j. In this context this means that the splitting procedures converge also 
in the operator norm and we get a uniform error estimate. It is a very interesting and 
promising direction, but unfortunately not applicable for a large class of equations, see the 
delay equation example in the last section. This is why we investigate strong convergence 
in this paper. 


3. Convergence of the splitting procedures 


First we show that the investigated splitting procedures are convergent if the exact 
solutions of the sub-problems are known. Then the convergence is showed also for the 
case when the solutions of the sub-problems are only approximated by certain spatial and 
temporal numerical schemes. 


3.1. Convergence with exact solutions. Since operator splitting procedures are finite 
difference methods, their convergence is defined as in Definition 12.31 More precisely, a 
splitting procedure is called convergent at a fixed time level t > 0 if 

lim || [F spl (t/n)] n x — U{t)x\\ = 0 for all x G A", 

n—>oo 

with index ‘sq’, ‘St’, or ‘w’ for the sequential, Strang, and weighted splitting, respec¬ 
tively. The convergence of the splittings will be a consequence of Chernoff’s Theorem (see 
Chernoff [(5}, and Engel and Nagel (j^, Corollary III.5.3]). Let us observe that Chernoff’s 
Theorem and one direction of Lax’s Theorem 12.61 state the same: the stable and consistent 
methods are convergent. 

Since the stability condition appearing in Chernoff’s Theorem will play an important role 
in obtaining the convergence of the splitting method, we cite Lemma 2.3 from the paper 


Lemma 3.1. Let us assume that there exist constants M > 1 and wet such that 
(11) || [S(t/n)T(t/n )] n || < Me ut for all t > 0, n G N. 

Then the followings hold. 

(i) There exist constants M\ > 1, oj\ G M such that 

|| [S(t/n)T(t/n)] n ~ 1 \\ < M x e Ult for all t > 0, n G N. 

(ii) There exist constants M 2 > 1, 1 x 2 G M such that 

|| \T{t/n)S(t/n)} 1 < M 2 e U2t for all t > 0, n G N. 

(in) There exist constants M 3 > 1, w 3 G R such that 

|| [S{t/2n)T{t/n)S{t/2n)] n \\ < M^e U3t for all t > 0, n G N. 

Lemma 13.11 states that the stability of the sequential splitting implies the stability of the 
Strang and of the reverse order splitting. 
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Remark 3.2. The stability assumption for the weighted splitting can be formulated as 
follows. There exist constants M 4 > 1, u 4 6 R such that 

(12) \\[eS(t/n)T(t/n) + {l-Q)T(t/n)S(t/n)] n \\ < M A e u * and 

|| \QT{t/n)S{t/n) + (1 - Q)S(t/n)T(t/n)\ n || < M 4 e“ 4t 

for all t > 0, n G N, where 0 G [0,1]. Notice that, in general Condition (ITTl) does not 
seem to imply the stability condition for the weighted splitting. The question of giving 
a simpler condition implying the stability condition for the weighted splitting is still an 
open problem. 

The following proposition is a consequence of Chernoff’s Theorem. 

Corollary 3.3. Under the Assumption 11.11 the sequential and the Strang splittings are 
convergent at a fixed time level t > 0 if the stability condition OT]) is satisfied. The 
weighted splitting is convergent if the stability condition (fT21) is satisfied. 

Proof. 

Sequential and Strang splittings: See [g, Corollaries 2.5 and 2.7]. 

Weighted splitting: 

In order to show the convergence, we apply Chernoff’s Theorem to the operator 
F w (h ) := QS(h)T(h) + (1 - Q)T(h)S{h) 0 G (0,1). 

The stability holds due to our assumptions. Using the semigroup property, we obtain that 

[QS{h)T(h) + (1 - Q)T(h)S(h)]x - x 
h —^0 h 

=0 ton mm*- x lim msw*- x 

h-> 0 h h-> 0 h 

=Q(Ax + Bx) + (1 — 0)(Ax + Bx) = (A + B)x 

for all x G D(A) D D(B) and 0 e (0,1), since the topology of pointwise convergence on a 
dense subset of X and the topology of uniform convergence on relatively compact subsets 
of X coincide on bounded subsets of C(X), see e.g., Engel and Nagel j§, Proposition A.3]. 
Thus, the consistency criterion of Chernoff’s Theorem holds as well. □ 

Remark 3.4. Notice that in the important special cases, where the semigroups S and 
T are quasi-contractions or commuting, all the stability conditions are automatically 
satisfied. 

3.2. Convergence with approximations. In the previous section we investigated the 
convergence of the splitting procedures in the case when the sub-problems are solved 
exactly. I 11 concrete problems, however, the exact solutions are not known. Therefore 
the use of a certain approximation scheme is needed to solve the sub-problems. When a 
partial differential equation is to be solved applying a splitting procedure together with 
approximation schemes, we have to follow these steps. 

1. The spatial differential operator is split into sub-operators of simpler form. (For 
instance according to the different physical phenomena or space directions, etc.) 

2. Each sub-operator is approximated by an appropriate spatial discretization scheme 
(called semi-discretization). Then we obtain systems of ordinary differential equa¬ 
tions corresponding to the sub-operators. 

3. Each solution of the semi-discretized system is obtained by using a time-discreti¬ 
zation method. 
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In this section we investigate the case when the semi-discretized systems can be solved 
analytically, i.e., without using a time-discretization method. That is, we assume that 
the semigroups are approximated by approximate semigroups (step 1 and step 2). In the 
second case the semigroups are approximated by operators which are not necessarily semi¬ 
groups. They represent the case when the solutions of the semi-discretized systems are 
obtained by using a time-discretization scheme (step 1, step 2, and step 3). This case will 
be investigated in our forthcoming work. We remark that the convergence of the splitting 
together with the time-discretization method (without the spatial approximation scheme: 
step 1 and step 3) is investigated, e.g., in [7]. 


Consider the abstract Cauchy problem (lACPj) on the Banach space A" for the sum of 
the generators (A, D(A)) and [B, D(B)y As introduced e.g. by Ito and Kappel in [TEI, 
Section 4.1] and Pazy in [H, Section 3.6], we define approximate spaces (“mesh”) and 
projection-like operators between the approximate spaces and the original space X. 

Definition 3.5. Let X m , m e N be Banach spaces and take operators 

Pm • t X rn and J m . X rn ^ X 

fulfilling the following properties: 

(i) P m J m = I m for all m e fd, where I m is the identity operator in X mi 

(ii) lim J m P m x = x for all x e X, 

m—>• oo 

(iii) || J m || < Mj and ||P m || < Mp for all m 6 N and some given constants Mj, Mp > 0. 


The operators P m together with the spaces X m usually refer to a kind of spatial dis¬ 
cretization (triangulation, Galerkin approximation, Fourier coefficients, etc.), the spaces 
X m are usually finite dimensional spaces and the operators J m refer to the interpolation 
method, how we associate specific elements of the function space to the elements of the 
approximating spaces (linear/polynomial/spline interpolation, etc.). 

First we split the operator A + B appearing in the original problem (lACPjl into the sub¬ 
operators A and B. In order to obtain the semi-discretized systems, the sub-operators 
A and B in equation (I ACPI) have to be approximated by operators A m and B m for 
m e N fixed. Assume that the operators A m and B rn generate the strongly continuous 
semigroups (T m (i)) and ( S m (t )) t>0 on the space X m , respectively. For the analysis of 
the convergence, we need the following definitions. 


Definition 3.6. We consider the following properties of the semigroups (T m (f)) t>0 , 
{S m (t)) t>Q , me N, and their generators (A m ,D(A m )), (B m ,D(B m )), me N, respec¬ 
tively. 

(i) Stability: 

there exist constants M, Mp, Ms > 1 and oj,ojt,ujs £ ® such that 

(a) ||T(h)|| < M T e UTh and ||T m (h)|| < M T e“ Th , 

(b) \\S(h)\\ < M s e^ h and \\S m {h)\\ < M s e^ h , 

for all h > 0, and either 

(13) \\[S m (t/n)T m (t/n)] k \\ < Me kul ™ for all t > 0, k,n,m e N 

in case of the sequential and Strang splittings, or 

[eS m (t/n)T m (t/n) + (1 - e)T m (t/n)S m (t/n)] k || < Me k ^ 


(14) 
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for a 0 G [0,1] and for all t > 0, k,n, m G N in case of the weighted splitting. 

(ii) Consistency: 

(a) lim J m A m P m x = Ax for all x G D(A), 

ra—>• oo 

(b) lim J m B m P m x = Bx for all x G D(B). 

m—too 

The semigroups T m , S m , m G N, are called approximate semigroups , and their generators 
[A m , D^Am )), (B m ,D{B m )) are called approximate generators if they possess the above 
properties. 

Corollary 3.7. From the assumptions in Defjnition \3.b\ and from the Trotter- Kato Ap¬ 
proximation Theorem (see Ito and Kappel lla . Theorem 2.1]) it follows that the approxi¬ 
mate semigroups converge to the original semigroups uniformly, that is: 

Convergence: 

(a) lim J m T m (h)P m x = T{h)x Vx G X and uniformly for h G [0, £ 0 ]> 

m—>• oo 

(b) lim J m S m (h)P m x = S{h)x Vx G X and uniformly for h G [0, £ 0 ]. 

m—>• oo 


We remark that the stability condition (i)(a) or (b) in Definition 13.61 is fulfilled for many 
cases in applications, for example, for the semidiscrete Galerkin method for parabolic 
equations (see Larsson and Thornee [13, Section 10.1]) or a large class of Galerkin t ype 
approximations for hyperbolic or more general equations (as an example, see Fabiano [lOj] 
or Fabiano and Turi |ll|). 


Definition 3.8. For the case of spatial approximation, we define the split solutions of 


(BCPD 


as 


< p Lm ■= j m K p, (t/n)rp m x 

for m, n G N fixed and for x G X, where index ‘sph is ‘sq’, ‘St’, or ‘w’. The operators 
F m , describing the approximation schemes together with the splitting procedures have 
the form 


(15) F^(h) := S m (h)T m {h), 

(16) F%(h) := T m (h/2)S m (h)T m (h/2 ), 

(17) FZ(h) := eS m (h)T m (h) + (1 - Q)T m (h)S m (h), 0 G (0,1) 

for the sequential, Strang, and weighted splittings, respectively, with h G (0, t 0 ] • 

Definition 3.9. The numerical method for solving problem (1ACPD described above is 
convergent at a fixed time level t > 0 if for all £ > 0 there exists IV G N such that for all 
n,m > N we have 

||< P L(X> ~u{t)\\ < £, 

where the index ‘sph refers to ‘sq’, ‘St’, or ‘w’, respectively. This is the usual well-known 
notion of the convergence of a sequence with two indices and we will use the notation 

lim = u{t) 

n,m—>oo ’ 

to express this. 

Remark 3.10. Observe that, analogously to the case of exact solutions (i.e., splitting 
without approximation), the stability condition (TT3l) implies the stability of the reversed 
order and the Strang splitting. Since the proof is analogous to the one of Lemma 13.11 we 
omit it. 
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This means that the sequential and Strang splittings fulfill their stability condition if the 
stability condition fTEill of the sequential splitting (with approximations) holds. Therefore, 
in this case, it suffices to control only this condition. Notice, however, that for the weighted 
splitting we need the more complicated condition (fTTl) (with approximation). 

Notice that the convergence of a sequence of two indices is in general difficult to treat, 
except in a well-known special case we recall here in a form we shall need in our proofs 
later on. 


Lemma 3.11. Let X be a Banach space, D C X a dense subset in X, and V(h), V m {h ) : 
D —> X operator for h G [0, to] an d V : D —> X. Let us assume that 

(18) 3 lim V m (h)x = V{h)x forallxED, uniformly for h G [0, to] and 

m—>oo 

(19) 3 lim V ( h)x = Vx for all x G D. 

h—>0 

Then 

lim V m [t/n)x = Vx holds for all x G D 

n,m—> oo 

and for all t G [0, to] fixed. 

In order to prove the convergence of operator splitting in this case, we state a modified 
version of Chernoff’s Theorem, which is applicable for approximate semigroups as well. 
The version we present here is a slight modification of Pazy j23|, Theorem. 6.7], which is 
modified to fit our setting. 


Theorem 3.12 (Modified Chernoff’s Theorem). Consider a sequence of functions F rn : 
1R + £( X m ), m G N, satisfying 

(20) F m ( 0) = I m for all m G N, 


and that there exists a constant M > 1 such that 

(21) HKWlInpM < M for all t > 0, m, k G N. 

Assume further that 


3 lim 

m—>oo 


uniformly in n GN and that 


( 22 ) 


Gx := lim lim 

n—>oo m— KX) 


Jm,Fm(t/rf) P m X J rn P rn X 
t/n 

JmFmif / nf P m X J m P m X 


t/n 


exists for all x E D C X and an arbitrary t > 0, where D and (A 0 — G)D are dense sub¬ 
spaces in X for Ao > 0. Then the closure G of G generates a bounded strongly continuous 
semigroup (U(t)) t>0 , which is given by 


(23) 


U{t)x — lim J m {F m (t/n)] n P m x 


n,m—>oo 

for all x G X uniformly for t in compact intervals. 


Proof. For h > 0 define 


G 


h,n,m 


F m (h/n ) 


h/n 


e C(X r , 
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for all n, m E N. Observe that lim J m Gh,n,mPmX = Gx for all x E D. Every semigroup 

( e t t>Q 


JG h , n ,m)^ n satisfies 


(24) 


?tGh ,1 


t k n k 


< e~ tnlh \\e^ F m(h/n)/h\\ < V' ——1| [F m (h/n)] fc || < M 

II II /. —/ h k k\ 


k =0 


for every t > 0. This shows that the assumptions of the Trotter -Kato Theorem (see Ito 
and Kappel (THI . Theorem. 4.2] together with Engel and Nagel jl, Theorem 111.4.8]) are 
fulfilled, and we can apply it first taking limit in m —> oo (which is uniform in n E N), and 
then taking limit as n —» oo. Hence, the closure G of G generates a strongly continuous 
semigroup U given by 

(25) lim \\U(t)x — J m e tGt ’ n ' m P m x || = 0 for all x E X 

n,m—>oo 


uniformly for t in compact intervals. On the other hand, we have by Lemma 111.5.1. in 
Engel and Nagel 


(26) 


| \J m e‘ G, '"P m x - J m \F m {t/n)}''P m x || 

-111 p n(Fm(t/n)-Im) p T _ T \ p (f/r>)] n P 

11 mX' 1 m\ L / 1 \ 1 n 

Jm.F'mif'/Tl)PmX ' 1 rnP'm,^\\ 


,X\ 


tM 


Jm,Pm(t/Tl) P m X 


n,m—>oo 


» 0 


t/n 

for all x E D, and uniformly for t in compact intervals. The combination of and 
yields 


\\U(t)x-J m [F m {t/n)] n P m x\\ 

<|| U(t)x - J m e tGt ’ n ' m P m x\\ + || - Jm[F m (t/n)\ 


PmX\\ -1 0 


for all x E D, and uniformly for t in compact intervals. By the uniform boundedness 
principle the statement follows for all x E X. □ 


Remark 3.13. Theorem 13. 1 21 remains valid in the case when the stability conditions 
reads as 

(27) \\[F m (t)] k \\ < Me kut 

for all t > 0 and m E N, k E N, and for some constants M > 1, w 6 1, as we can see it 
from a standard rescaling procedure, see Engel and Nagel jl^Corollary 111.5.3]. 

In the remainder of this subsection, we consider the convergence of the various splitting 
procedures. 


Lemma 3.14. Let J m , P rn , T m be operators defined in Definitions 1.7,51 and \S. 61 Then the 
following holds for all t E [0, t 0 ]-' 

JrnPmifi / Tl) PmX JmPmX 


lim 

n,ra—>-oo 


t/n 


= Ax for all x E D(A ), 


where the limit as m -E oo is uniform in n gN. 

Proof. It suffices to prove that the operators 

dmPm ( h ) Pm JmPrr 


V m {h) = 


h 
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fulfill the conditions in Lcmma l3.Hl In order to determine the first limit f[T8j) . the following 
norm has to be investigated for all x G D(A): 


Jmijl) Pm-E ^rnPnnA T^PjX X 


h h 

I J m A m T m (s)P m xds - / HT(s)xds 


h h 

o o 

^ SUp || JmA m T m (s') P m X AT(s):r|| SUp || JrnPm T(s)A:r|| 

sG[0,to] sS[0,lo] 

< sup II J m T m (s)P m [J m A rn P m x - Ax] + [J m T m (s)P m - T(s)]Ac|| 

se[o,to] 

< sup \\J rn \\\\T m (s)\\\\P m \\\\J m A m P m x - Ax\\+ sup \\[J m T m (s)P m -T(s)]Ax\ 

se[0,to] sG [0,tol 


According to Definition 13.61 the term || J m A m P m x — Ax|| tends to 0 as m tends to infinity. 
Since y := Ax is a fixed element in the Banach space X, || J m T m (s)P m y — T(s)i/|| tends 
to 0 uniformly in h as m tends to infinity because of Corollary 13.71 Operators J m and 
P m were assumed to be bounded. Term ||T m (s)|| was also assumed to be exponentially 
bounded with a constant ojt, which is independent of m. Therefore, 

sup ||T m (s)|| < sup M t q U tS < Mre max{0,a,T} *° = const. < oo. 

sG[0,to] s€[0,to] 


Hence, the above difference tends to 0 uniformly in h, therefore, (1T5|) holds. The second 
limit (ITU]) can be obtained by using the definition of the generator: 

lim - ■}_ _ Ax y x e 

o h 

Thus, the statement is proved. □ 


Corollary 3.15. The following statement can be proved similarly as Lemma \3. If 

/ AjP m X J m P rn X 


lim 

n.m—>oo 


t/n 


= Bx 


Vx e D(B). 


Now we really turn our attention to the convergence of the different splitting procedures. 
First we show the convergence of the split solution defined in (1 151) . i.e., in the case when 
the sequential splitting is applied. 

Theorem 3.16. The sequential splitting KT5 1) is convergent at time level t > 0 if the 
stability condition holds for the approximate semigroups, and the approximate generators 
are consistent according to Definition ^. (X 


Proof. According to Chernoff’s Theorem 13. 121 the sequential splitting is convergent if the 
stability (TUT]) and the consistency (122]) hold for the operator 


(28) F m (h) := S m {h)T m {h). 

The stability condition (1271) is fulfilled, since we assumed that (fT3j) holds. In order to prove 
the consistency criterion ([22]) . the following equation should be verified: 


lim 

n,m—>oo 


JmBmif / Tl) P'm.'E JmPm 


X 


t/n 


(.A + B)x for all x e D(A) D D(B). 
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We investigate the following limit: 


lim 

n.m— 


J m*5 m(t / Af)T m (t/Tl) P m X Jm/PmA 


lilll < JmSm^t / Tl) P r 

n,m—KX) 


t/n 

JmTm(t/n)P m X JmPm% 


t/n 


lim 

n.m—KX5 


Jm,Sm(t/’fl) Pm% JmPm% 
t/n 


where we apply Lemma 13.111 for each term. Choosing 


Vm(h) ■ JmPmiJPjPm 
and applying Corollary 13.71 we obtain that 

lim J m S m (h)P m x = S(h)x Vx G X uniformly for h G [0, to] and 

m—>■ oo 

lim S(h)x = x for all x G X. 
h->o 

Notice further that the set {^(J m T m (h)P m x — J m P m x ) : h G [0, to]} is relative compact 
for x G D(A), and that on compact sets the strong and the uniform convergence is 
equivalent. Hence, Lemma 13. 141 and Corollary 13.151 can be applied, thus, we obtain that 


(29) 


lim 

n,m—>oo 


= lim 

n.m—KX) 


i/ m P m (t/?7.)P m X J m P m X 
t/n 

JmSmi/// fl) PmX J m P m X 
t/n 


{I A + B)x 


for all x G D(A) D D(B). Here / G B(X) denotes the identity operator. This means that 
the consistency criterion of Chernoff’s Theorem is fulfilled for the sequential splitting, 
thus, it is convergent with spatial approximations. □ 


Now we prove the convergence of the split solution defined in (fTBlh i.e., in the case when 
the Strang splitting is applied. 

Theorem 3.17. The Strang splitting (ITU is convergent at time level t > 0 if the stabil¬ 
ity condition holds for the approximate semigroups, and the approximate generators are 
consistent according to Definitio ns. (A 

Proof. As shown in Lemma 13.11 the stability condition for Strang splitting follows from 
that of sequential splitting which has already been assumed in Def ini tion 13.61 Hence, now 
it suffices to check whether the consistency (l22|i holds for the operator 


F m (h) := T m (h/2)S m (h)T m (h/2). 
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Thus, we should investigate the following limit 

Jm,Fm(f / F) P m X Jm,Fm% 


lim 

n,m—> oo 

= lim 

n,m—> oo 


t/n 

J m T m (t/2n)S m (t/n)T m (t/2n)P m x - J m PmX 


t/n 


= lim J m T m (t/2n)S m {t/ri)P 1 

n,m—>oo 


JmTm{t / '2) P m X J m P m X 


lim J m T m (t/2n)P 7 


t/n 

dm{t /' 2 ) PmX JmPm% 


1 1 TT1 J P 
11111 'J m. 1 rr 


t/n 

dm,Fm(t / '2) PmX J m P m X 
t/n 


We can apply Lemma [3.141 and Corollary 13.151 again to each term, and obtain that 


lim 

n,m—Kx) 

= lim 

n,m—>-oo 


Jm,Pm{t / Tl) P m X J m P m X 
t/n 

JmT m (t/2n)S m (t/n)T m (t/2n)P m x - J m P m x 
t/n 


=/ ■ / ■ \Ax + / • Bx + iylx = (A + P)x 

for all x G D(A) D D{B). Therefore, the Strang splitting is convergent with spatial ap¬ 
proximations. □ 


Finally, let us prove the convergence of the split solution defined in (1171) . i.e., in the case 
when the weighted splitting is applied. 

Theorem 3.18. The weighted splitting [Tip is convergent at time level t > 0 if the stability 
condition m holds for the approximate semigroups, and the approximate generators are 
consistent according to Definition ^. (A 

Proof. Similar to the calculations presented in proofs of Theorems 13.161 and 13.171 for the 
operator 

F m {h) := eS m (h)T m (h) + (1 - Q)T m (h)S m {h), 0 G (0,1). 

□ 


Summarizing, we can say that by Chernoff’s Theorem 13. 121 the stability condition (1IB|) or 
(ED implies the convergence of the splitting procedures in the case of consistent approx¬ 
imations. 


4. Operator splittings for delay equations 

Since many physical processes depend on a former state of the system as well, they have 
to be described by delay differential equations containing a term depending on the history 
function (see Q, Kappel [l3|). Although these differential equations cannot be written as 
an abstract Cauchy problem on the original state space X, their solutions can be obtained 
by an operator semigroup on an appropriate function space (called phase space). For a 
systematic treatment of the problem we refer to the monograph [4j. We concentrate here 
on the case of bounded delay operators, while unbounded delay operators will be treated 
in a subsequent work. 
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Consider the abstract delay equation in the following form (see, e. 


g; fl): 


(DE) 


= Cu(t) + t> 0, 

w(0) = x E X, 
u 0 = f EL 1 (1-1 : 0},X) 


on the Banach space A", where (C, D{C)) is a generator of a strongly continuous semi¬ 
group on X, and $ : L 1 ([—1, 0], AC) —* X is a bounded and linear operator. The history 
function u t is defined by u t (a ) := u{t + a) for a E [—1, 0]. 


In order to rewrite (IDE!) as an abstract Cauchy problem, we take the product space 
£ := X x L 1 ([— 1, 0], A") and the new unknown function as 


1 1— y C/(f) •— 


u(t) 

Ut 


E £. 


Then (IDED can be written as an abstract Cauchy problem on the space £ in the following 
way: 

d U(t) 


(MV) 


= GU(t), t > 0, 


dt 

U( 0) = (*) e E, 


where the operator (Q, D(Q)^) is given by the matrix 


(30) 

on the domain 


G ■= 


C $ 
0 -f 

da 


D(S ) := {© £ D(C) x wyi-l.o], X) : r,(0) = ?} . 


It is shown in [3, Corollary 3.5, Proposition 3.9] that the delay equation (IDED and the 
abstract Cauchy problem <\ACV\] are equivalent, i.e., they have the same solutions. More 
precisely, the first coordinate of the solution of (I ACV\) always solves (IDEE Due to this 
equivalence, the delay equation is well-posed if and only if the operator (G, D(G)) gener¬ 
ates a strongly continuous semigroup on the space £. 


The following case was partly investigated in the paper 

Assumptions 4.1. 

(a) The operator (C,D(C)) generates a strongly continuous contraction semigroup 
TW) t >„ 0.1 A'. 

(b) The delay operator <h : L 1 ([—1,0], A) —y X is bounded. 

Since the delay operator $ is bounded, the delay equation (IDEI) is well-posed by 
Theorem. 3.26]. In order to apply an operator splitting procedure, we split the operator 
in (I ACPI) as 


G = A + B, 
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where the sub-operators have the forms 


(31) 


A:= 
B : = 


C 0 \ 

o £) 
0 $ \ 

0 0 ) ’ 


D(A) := D(G), 
D(B) :=£. 


Since C is a generator and $ is bounded, the operators A and B generate 
continuous semigroups (T(i)) (>0 and (<S(i)) , respectively. It is shown in 


3.25] that T is given by 


( V(t) 

V V t T 0 (t) 


1 


the strongly 
||_4, Theorem 


where (T 0 (i)) t>0 is the left shift semigroup defined by 




0, 


for all / G L p ([—1, 0], A") , and V t is 


{V t x){<r) 


V(t + a)x, 
0, 


if a G [—1, —t), 
if a G [—t, 0], 


if a G [—t, 0], 
if a G [—1, —t). 


for all x G X. Since <3> is a bounded operator, B is also bounded on £. Therefore, the 
semigroup S generated by B is 


S(t) := e tB = 1 + tB 


1 M> \ 

0 I ) ’ 


where /, /, and 1 denote the identity operators on X , L 1 ([—1, 0], X), and £, respectively. 
By formulae (|6|). (151), and (TTO]) of the sequential, Strang, and weighted splittings, the split 
solutions of the delay equation with initial value (J) G £ p can be written as 

(32) tC(«) = [5(«/n)T(t/r.)]”(J), 

(33) Ul\t) = [T(t/2n)S(t/n)T (t/2ji)] n (*), 

(34) M-(t) = [eS(t/n)T(t/n) + (1 - e)T(t/™)S(t/r,)]”$ 

for n G N hxed and 0 G (0,1), for the sequential, Strang, and weighted splittings, 
respectively. 

Theorem 4.2 (Theorem. 4.2, Corollary 4.3 in (<§]). Under the Assumptions \4-l\ the se¬ 
quential, Strang, and weighted splittings applied to the delay equation urn with sub¬ 
operators (EIP are convergent at a fixed time level t > 0. 


By Proposition 13.31 we only have to show that the stability condition (HIT) is fulfilled. For 
the proof, see [8j. (We note that the stability of the weighted splitting follows from the 
estimates in |8| because we can choose M — 1 for both semigroups, see Remark 13.41 1 
In Theorem 14.21 we showed that the sequential, Strang, and weighted splitting procedures 
are convergent when they are applied to the abstract Cauchy problem (L4C'PI) associated 
to the delay equation (1DED . We now combine these results with spatial approximations. 
Similarly as in Section 13.21 we define the following spaces and operators. 


Definition 4.3. For m G N we take 
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(i) X m Banach spaces, 

(ii) P m and J m operators satisfying the conditions in Definition 13.51 for the Banach 
spaces X and X m , 

(iii) P m and J m operators between the Banach function spaces L 1 ([—1,0], X) and 
L 1 ([— 1, 0], X m ) defined as 

( Pmf){<7 ) := P m f(cr) and (J m / m )(a) := J m fm(cr) 

for all / G L 1 ([—1, 0], X) and f m G ^([-1, 0], X m ), 

(iv) the spaces £ X x L 1 ([—1,0],X) and £ m := X x L 1 ([— 1, 0], X m ) , and 

(v) the operators 

Vm := ( P o L ) iu £ ^ Jm - = ( J o" £) iu 


Proposition 4.4. The operators V mi J m satisfy the approximating properties defined in 
Definition 15.51 


Proof. • \\V m \\ < Mj, since: 


I IP™ 11 = max 
IWI<LII/I 


<i 


^x\\,\\ p mf\\} < „ max ] 


sup \\P m f(s)\\} 
o-e[—i,o] 


< max {||P m ||||x||,||P m || sup ||/(s)||} < max{||P m ||, ||P m ||} 
M |< 1,||/||<1 o - g [— 1 , 0 ] 

= ||Pm|| < Mp, 


and 11 J m 11 < Mj similarly. 

• PmJm = P-m, where T m denotes the identity operator on £ m , since 


f Pm ® \ ^ / PrrJrn _ \ 

\ 0 Pm ) V 0 Jrn ) \ 0 P m J m J 


where P m J m = I m and (P m J m / m )(s) 

^([-1,0],X m ). 


• lim J m V m ( x f ) = (f), because 

nn—too 7 


PmJmfm ( & ) 


fm($fi hence PmJm fn in 


J m Pm% X 

JmPmf - / 


= max { \\J m P m x - x||, II JmPmf - / lllA 


where J m P m x converges to x because of Definition 13.51 and 


\\JmPmf-f\\u = 

0 0 

= j | \(JmPmf ~ /)(cr) ||dcr = j || J m Pmf (cr) ~ /(<t) 11 dcr ———> 0, 


-1 


-1 
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by Lebesgue’s Theorem since it converges to zero for all a G [—1, 0] and J m Pmf — / 
is bounded: 

o 

\\JmPmf ~ f\W < \\JmPmf\\v + II/IIl 1 < J || (<?)\\d(T + \\f || L i 

-1 

0 0 
= j \\ JmPmf {o-)\\da + ||/|| L 1 < II J m ||||P m || J Wf(a)\\da+\\f\\ L i 

-l -l 

<MjMp\\f\W + II/IIli = ( MjMp + 1)||/|| L 1. 

□ 

From now on we assume the following, using the notations and terminology of the 
previous section (see Assumptions m- 

Assumptions 4.5. 

(a) There exist operators C m : X m —> X rn generating strongly continuous semigroups 
V m satisfying ||V^ n (t)|| < 1, such that J m C m P m x —> x for all x G A". 

(b) There are bounded operators <3> m : L 1 ([—1, 0], X m ) —> X m satisfying ||$ m || < c||$||, 
such that Jm^mPmf $/ for every / G L 1 ([-1, 0], X). 

Remark 4.6. A typical and usual choice for the operators <f> m would be the one obtained 
through the spatial discretization as $ m / m := P m ^Jmfm- 

We are now in the position to define the approximating operators for the delay equation. 
Consider the operators 


C m 0 


da 


(35) Am ■■= 

on the domain 

D(Am) := {(y G D(C m ) x Wfo([—1,0], X m ) : r, m ( 0) = , 


0 <f> r 


and the operators 

y J , D(B m ) := S m , 

corresponding to the spatially discretized (ordinary) delay equations 

< dw m (t) 


(DE r 


dt = C m u m (t ) + $ m u?, t > 0, 

U m ( 0 ) = x m G X m , 


U r 


= / m GL 1 ([-l,0],A m ). 


Theorem 4.7. Under these assumptions, the sequential, the Strang, and the weighted 
splittings are convergent for the delay semigroup. 


Proof. By Theorems 13.16113.171 and 13.181 we have to show that the stability and consis¬ 
tency conditions of Definition 13.61 are satisfied for the operators A m , B m . 
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It is clear that the operators A m generate strongly continuous semigroups 


T m (t) := 


Vjt) 0 \ 

V t m T 0 (t ) ) 


where 


(V t m x){o) 


V m (t + a)x, 

0, 


if cr G [—t, 0], 
if cr G [—1, —t). 


From the calculations in the proof of [N. Theorem 4.21, especially formulae 181/(15) and 
0/(16), we deduce 


T m (f)\\ < 1 + t, 


and further, by the assumptions on the operators Q m , B m generate norm-continuous semi¬ 
groups 



satisfying 


\\S m (t)\\ < 1 + tc||$||. 


Hence, the stability conditions (i)(a) and (b) from Definition 13.61 are satisfied. Further, 
from the same arguments as in the proof of [8j, Theorem 4.2] we obtain that the stability 
conditions (1T5D and fflTD are satisfied with M — 1 and uu — 1 + c|| < F||. Here we also have to 
use the fact that t/n denotes the time step, and k is the number of steps we have already 
done in order to reach the time level t. This means that the relation k < n always holds. 
Then we can apply the following estimate: 


\\[S m {t/n)T m {t/n)] k \\ < \\[S m (t/n)T m (t/n)] n \\ < \\S m {t/n)\\ n \\T m (t/n)\\ n < e^^"/ 


Now we only have to check the consistency conditions. The condition (ii)(b) in Definition 
13.61 follows immediately form our Assumption 14.51 (b). 

For condition (ii)(a), observe first that if (j) G D(A), then /(0) = x, implying P m f( 0) = 
P, m x , and so (p m p G D(A m ). Hence, V rn D(A) C D(A m )- Further, for (*) G D(A), we 
have that 

'TAT) ( x \ T A (P nX \ s~r (Crn,Pm%\ / JrnCrnPrn^\ 

Jm^rnt m{ f ) ~ Jm^m\ Prn f) ~ “ V )' 

Again, by our assumption on the operators C m and by Lebesgue’s Theorem we see that 

7 A V (D — 5 . ( Cx ) 

sj rri' AX m > m \ jj 7 \ f' J ‘ 


Hence, the desired stability and consistency conditions are satisfied meaning that the 
splittings are convergent. □ 
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